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Abstract 

Thermalization of a heavy quark near rest is controned by the correlator of two electric fields 
along a temporal Wilson line. We address this correlator within real-time, classical lattice 
Yang-Mills theory, and elaborate on the analogies that exist with the dynamics of hot QCD. 
In the weak-coupling limit, it can be shown analytically that the dynamics on the two sides 
are closely related to each other. For intermediate couplings, we carry out non-perturbative 
simulations within the classical theory, showing that the leading term in the weak-coupling 
expansion significantly underestimates the heavy quark thermalization rate. Our analytic 
and numerical results also yield a general understanding concerning the overall shape of 
the spectral function corresponding to the electric field correlator, which may be helpful in 
subsequent efforts to reconstruct it from Euclidean lattice Monte Carlo simulations. 
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1. Introduction 



The real-time dynamics of heavy ion colhsions is governed by QCD at relatively large cou- 
pling, which remains poorly understood despite significant theoretical efforts. While in prin- 
ciple both the weak-coupling expansion and lattice QCD provide systematically improvable 
schemes for the calculation of any physical quantity, including unequal-time correlation func- 
tions (At ^ ^/T) at a finite temperature T (T > 200 MeV), both are in practice faced with 
serious limitations, related to the convergence of the weak-coupling expansion and to the need 
to carry out analytic continuation, respectively. Therefore many current attempts to describe 
the real-time dynamics of QCD at realistic temperatures rely on models. For instance, for 
heavy quark thermalization and diffusion, the main topics of the present paper, a relatively 
successful model treatment can be obtained by incorporating bound states as dynamical de- 
grees of freedom into the description of an otherwise partonic medium [Ij. Unfortunately, 
such models need typically to be tuned to the particular observable in question, rather than 
having a universal character, and also do not allow for a systematic improvement. 

As an alternative to models, much recent literature has focused on analogue theories, 
by which we mean well-defined frameworks which are sufficiently close to QCD that most 
interesting QCD measurables have equivalents in the analogue, but which are nevertheless 
more amenable to calculation. Two such frameworks have been especially widely used: QCD 
truncated to the first non-trivial order in the weak-coupling expansion (see, e.g., ref. [2]), and 
A/'=4 Super- Yang- Mills theory in the limit of an infinite number of colors and a large 't Hooft 
coupling (see, e.g., refs. [3l|4j). 

QCD truncated to the first non-trivial order in the weak-coupling series is the starting 
point of a systematic expansion, and thus arguably the most similar analogue theory, guar- 
anteed to be correct in the limit of a high temperature. The problem is that in the present 
setting it is technically extremely hard to work out subsequent terms in the weak-coupling 
series, given that extensive resummations are needed for dynamical quantities evolving over 
long time scales, and that in general several terms in the expansion would be needed, in 
order to obtain any kind of convergence. In fact, even though five subleading orders are 
available for thermodynamic (i.e. equal-time) quantities such as the pressure [5], convergence 
remains debatable [6]; in the dynamical case at most the first non-trivial order has been 
reached so far [3 El [9l [TO] , and the results certainly display very large 0{g) corrections. This 
appears to indicate that the weak-coupling expansion is not well behaved except at very high 
temperatures, casting doubts on the physical relevance of truncated results in the realistic 
temperature range. 

Super- Yang-Mills theory (SYM) resembles QCD in that it is a gauge theory with matter. 
However it contains many more matter multiplets than ordinary QCD, and in a different 
representation of the gauge group, which makes the matching between the theories ambigu- 
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ous [TT]. Furthermore, computations on the SYM side are simple (using the famous AdS/CFT 
correspondence |12j ) only in the limit of an infinite coupling, whereas the interesting regime 
is probably intermediate coupling. 

In this paper we argue in favor of another analogue theory for the real-time dynamics of 
QCD: classical Yang-Mills theory regulated on a spatial lattice. This theory was developed 
by Kogut and Susskind [T3] and has been used extensively to study the rate of Chern-Simons 
number diffusion in Yang- Mills theory [H] , as well as (partly) out-of-equilibrium phenomena 
such as plasma instabilities [15], the dynamics of electroweak symmetry breaking [16], and 
inflationary preheating [17]. Recently, it was also applied for estimating the imaginary part 
of the real-time heavy-quark potential in QCD [18], and analogous methods were used for 
studying jet energy loss and transverse broadening in a hot non-Abelian plasma [19j. In this 
approach the infrared (IR) behavior of QCD is approximated by introducing semiclassical 
fields, while in the ultraviolet (UV) the quantum mechanical "cutoff" on thermal effects from 
short distances is replaced with a lattice cutoff. Formally, the classical limit corresponds to 
taking ^ — > 0, which is a non-trivial limit at non-zero temperatures |201l21j. and non-singular 
in the presence of the lattice cutoff. 

In this paper we use this framework to study one of the simplest gauge-invariant observ- 
ables, the correlator of two electric fields along a Wilson line, 

5 Ei=i / cite*"* Tr (?7(-oo - t) gEi{t, 0) U{t, 0) gEi{^, 0) C/(0, -oo)) 

"^'^^ " T¥ (C/(-oo - i/?, t)U{t, 0)C/(0, -oo)) ' ^^'^^ 

where (3 = 1/T; gEi = i[DQ,Di] is the color-electric field; U{ti),ta) represents a temporal 
Wilson line from ta to t{, at a fixed spatial location x = 0; and the trace is over the funda- 
mental representation. The denominator removes any regularization issues associated with 
the Wilson lines themselves. The zero-frequency limit «;(0) of this correlator is the momen- 
tum diffusion coefficient of a heavy quark [U [22], and the combination r]£i = K{0)/2M]^iy^T 
emerging from linear response relations, with Mkin denoting the so-called kinetic mass of the 
heavy quark, determines the heavy quark thermalization rate [2t I22|, [23]. 

Our goals and the organization of the paper are as follows. In Sec. [21 we describe the basic 
ideas behind classical lattice gauge theory as a tool for studying real-time quantities in QCD. 
In Sec. [3l we focus more precisely on the observable in Eq. (jl.ip . and use the weak-coupling 
regime to study how close the analogy between the two theories really is. The limitations 
of classical lattice gauge theory as a model for QCD are also illustrated, by studying the 
unphysical "strong-coupling" limit of the lattice-regulated theory. In Sec. [U we discuss 
the results we obtain at intermediate couplings, where the weak-coupling expansion fails 
yet classical lattice gauge theory still captures the correct infrared dynamics that causes the 
failure. Some discussion and conclusions can be found in Sec. [5l while two appendices contain 
details related to the weak and strong-coupling regime on the classical lattice, respectively. 
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2. Classical Lattice Theory: Basic Idea 

At a temperature far above the confinement scale, such that the effective gauge coupUng 
g is small, QCD (whether pure-glue or with dynamical quarks) possesses three different 
parametric length scales (= inverse momentum scales): 

• the length scale (vrT)"^, where most of the energy resides; 

• the "color-electric" length scale {gT)~^, where plasma screening effects become impor- 
tant and perturbation theory needs to be resummed [24|; and 

• the "color-magnetic" length scale (g^T/vr)^^, where interactions become genuinely non- 
perturbative. The longest spatial correlation lengths of gauge invariant operators are 
on this scale [25] . 

Different real-time correlation functions and physical properties of QCD are sensitive to 
different scales. For instance, the Chern-Simons diffusion rate is sensitive dominantly to the 
momentum scale g'^T/n; it depends on the other scales only in that they change the dynamics 
on this scale. Scattering, radiation and energy loss are sensitive mostly to the scale gT. This 
scale therefore captures much of the physics of current interest in heavy ion collisions, such 
as jet quenching and heavy quark thermalization. Shear viscosity, on the other hand, is 
principally sensitive to the scale ttT, since most of the energy and momentum reside there. 

Since for the momentum scales p ~ gT,g'^T/iT the loop expansion parameter related to 
bosonic fields, e ~ g'^h/{e^^P — 1), can parametrically be replaced with its classical limit, 
g'^T/p, it can be argued [26] and shown formally [20] that the physics at the scales gT, g'^T/ir 
is described by classical statistical field theory. Quantum mechanics is only relevant at the 
scale ttT, where its role is to ensure that thermal excitations on short scales are suppressed. 
The idea of the lattice analogue theory is to suppress thermal excitations on short scales 
instead by imposing a spatial lattice cutoff. The resulting theory is a classical field theory on 
all scales. At weak coupling it remains a three-scale theory, with: 

• the length scale a, where most of the energy resides; 

• the color-electric length scale {g'^T / a)^^/"^ ~ 0,(3^'^ , where plasma screening effects 
become important and perturbation theory needs to be resummed [271 128j: and 

• the length scale {g^T/ir)^^ ~ a/^L) where interactions become non-perturbative. The 
longest spatial correlation lengths of gauge invariant operators are on this scale. 

Here we have introduced the "lattice coupling", /3l = 2Nc/ g'^Ta, which controls whether 
interactions are perturbative at the lattice spacing scale. 

Physics at the scale o is definitely different from physics of the quantum theory at the scale 
(vrT)"^. In particular, lattice discretization breaks translational invariance so there is no 
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conserved momentum. This changes hydrodynamic behavior in an essential way, so one should 
not try to study shear viscosity with the classical theory. However, the more infrared scales, 
describing "collective phenomena", are only changed to the extent that the loop effects they 
feel from the hard momenta, associated with the expansion parameter e g^T/p g^Ta, 
differ from the corresponding quantum loop effects, with e ~ 5^/71". 

For equal-time quantities, these radiative effects turn out to be rotationally invariant and 
of exactly the same form as in the continuum quantum theory. They are simply a Debye 
mass parameter for the Aq field, of magnitude [29] 

"iD.cont = Q 9 T , (2.1) 

^DMt = ' S = 3.175911536... . 2.2 

47r a 

Equating these gives a concrete way of relating the lattice spacing a and the temperature T, 
a ~ 3TiNc/ttT{2Nc + A'^f). [Note that the p ^ 1/a lattice modes are playing the role both of 
ultraviolet gluonic degrees of freedom (the 2A'^c) as well as of quark fields (the A^f); infrared 
quarks can be neglected because at low frequencies the Fermi-Dirac distribution function is 
much smaller than the Bose-Einstein one.] At the dynamical level, however, radiative effects 
are no longer rotationally invariant [27\ I28j . which means that any "matching" between 
the lattice scale a and the temperature T is ambiguous, and only makes sense in order-of- 
magnitude. We return to this issue in more detail in the next section. 

Considering finally the color-magnetic scale, it remains the same, g'^T/n, in both theories. 
In other words, the precise form of the ultraviolet regulator is invisible to physics at the largest 
distances. At the same time, the dynamics on the scale g'^T/n is non-perturbative [25], but 
can relatively easily be simulated numerically through the classical description. 

Consider now the "weak-coupling regime" where the three scales are widely separated, and 
an observable dominantly determined by the scale gT. Its exact value is given by the leading 
order result modified by relative corrections suppressed by e ~ g'^T/p. While corrections 
from the hard scale p ~ ttT may remain controllably small down to low temperatures (cf., 
e.g., ref. [30j). experience with many observables such as the plasmon frequency [7|, the 
heavy quark diffusion coefficient [8j, the light quark dispersion relation [9], the jet quenching 
parameter q [10], or the Debye screening length [31j, has shown that radiative corrections 
from the scale p ~ gT itself can be very large (even if parametrically perturbative) . The 
challenge would therefore be to sum the corrections from the scales p ~ gT,g'^T/TT to all 
orders. 

The key property of the classical lattice theory is that it is amenable to a numerical 
simulation, and therefore indeed allows for all-orders resummations of the type mentioned 
to be carried out in practice. Should the large corrections come from the scale p ~ gT in 
the quantum theory, they are somewhat distorted in the classical lattice gauge theory, but 
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the results are still representative of the qualitative behavior. The contact to the quantum 
theory is only lost in the "strong-coupling regime" where Pi, <1; then all three scales are of 
the same order and the physics differs essentially from the quantum theory. 



3. Classical Lattice Theory: Electric Field Correlator 

In the previous section we argued on general grounds that classical lattice gauge theory and 
thermal QCD have qualitative and, in the very infrared, even quantitative similarities. We 
now want to demonstrate this explicitly for the case of heavy quark thermalization. 

Let us start by recalling the reason for why the zero-frequency limit of Eq. (jl.ip describes 
heavy quark momentum diffusion [4J and consequently thermalization [21 [22l [23] (on the 
formal level, the correspondence can be derived by making use of the heavy quark effective 
field theory [22]). In a classical framework, it is quite easy to see why this is the case. 
Intuitively, the field gEi exerts a Lorentz force on a charge carrier, and k = n{0) is the total 
correlation of that force with previous and future forces. The force changes the momentum 
of the heavy quarks. The classical lattice theory does not contain any heavy quarks, but 
we can still evaluate the force-force correlation function to see what momentum diffusion a 
heavy quark would feel. 

We discretize classical lattice gauge theory as in ref. fl4] and sample its thermal ensemble 
using the algorithm of ref. [32] (we recommend ref. [32] for a more detailed description of the 
procedure). The classical simulation is generally carried out in a gauge where the temporal 
links equal unity; thereby all the Wilson lines disappear from Eq. (jl.ip . and we only need to 
correlate the electric fields. For a comparison we have also carried out some simulations with 
the "improved" lattice action of ref. [33] . which provides a dispersion relation conforming 
more tightly to the continuum one, though it also differs significantly for p ~ 1/a. On the 
quantum theory side operator ordering plays a role; in the following we assume symmetric 
ordering (for details, see ref. [22]), whereby the quantum correlator has the same symmetries 
as the classical one; this can be obtained by k(uj) [^i^^] + '^(~'^)]/2 from the case literally 
shown in Eq. (II. ip . 

To start with we compare the theories at the free level. In the continuum theory, 

Kcont(w) = [1 + 2nB(a;)] ^ , (3.1) 

Dvr 



«^cont(i) = g^CFT^Tr"^ 



COSh^TTtT) 1 



(3.2) 



.sinh'^(7rtr) 3sinh^(7rtT)_ 
where is the Bose-Einstein distribution function. The vacuum behavior of /^-cont 

{t) reads 
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Figure 1: The correlation function K{t) (left) and its Fourier transform k{uj) (right) in the quantum 
continuum and classical lattice theories, at leading order (free level). To relate the theories, we convert 
the lattice spacing a to the inverse temperature 1 /T by equating Debye lengths for the pure-glue theory 
{N[ = 0), whereby a = 3S/27rr (cf. discussion below Eq. ^^). 



g'^CpT'^ir'^ / {irtT)'^ . In the classical lattice theory, on the other hand 



i^-utti^) = 7, / Kp-^)+Kp + ^) ; (3-3) 



2g^CFT rl'^ d3p 
'«iatt(t) = 7. / 7Tr^cos(pt), (3.4) 

where p = \/^, fP' = J2iPi = Yliia sin ^ )^. The results are plotteci^ in Fig. [TJ At first 
sight the lattice correlator and the vacuum-subtracted thermal correlator do not look alike: 
the main difference is in the large-time behavior, where the continuum correlator dies away 
but the lattice correlator displays decaying oscillations. The difference is explained when we 
look at the frequency-domain correlation functions. Here we see that the lattice correlator 
has cusps while the continuum correlator is smooth. The cusps are van Hove singularities 
which arise because the lattice excitations follow a modified dispersion relatiorl^, 

Wlitt(p)=p^ (3-5) 

which has vanishing slope at the corners of the Brillouin zone, p = {ni,n2,n^)'K / a, leading 
to cusps in the density of states at ojf^^^_ = (4, 8, 12)/a^. These van Hove singularities are 

^The "thermal part" is the difference of the full and vacuum parts, and it is this difference which is relevant 
for heavy quark thermalization. 

^The dispersion relation for the "improved" action is more complicated, see Eqs. (66,67) of ref. [33] 
[Eqs. (63,64) in the journal version]; the overall sign is wrong in the latter equation. 
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well understood and have little impact on the small-u; behavior which is actually of interest. 
In the small-frequency region to <^ T the theories do agree completely. In other words, the 
rather dramatic difference in the time-domain behavior shown on the left in Fig. [1] arises 
because n{t) is principally sensitive to the T or 1/a scale, where the theories are different; 
when we look at the frequency domain, the large frequency behaviors are very different as 
expected, but the low frequency parts agree. However, at the free level, the intercept k(0) is 
zero on both sides, so we need to consider interactions. 

The leading non-zero value for k(0) turns out to involve a logarithm of the scales vrT and 
gT. Considering the logarithm from the IR side, its origin lies in the fact that the electric 
gauge field self-energy gets an imaginary part for \uj\ < |p[ ~ gT, corresponding to the 
phenomenon of Landau damping (it also gets a real part, corresponding to Debye screening). 
The result for k is related to the cut of the electric field propagator, and can be written as 

S^g^CpNc f d^p p2 I- 



" 3 X« H2vr)3 ip2 + mlr J (2^ ^^^^ " '^^ " ^ ^ ^^^^^^^^ + "-(^^^ 

^'^-^<fln^ + ..), (3.6) 



67r \ rriD 

where Cp = {N^ — 1) /2A'c is the Casimir of the heavy quark representation. Here we omitted 
for notational simplicity quarks, and carried out the integral 2g'^Nc nB{q) [1+713(9)] = Trn^. 
The full computation and the result for the coefficient accompanying the logarithm can be 
found in ref. [2]. 

Consider then the classical lattice theory side. Restricting again to the leading logarithmic 
order, the only modifications needed are as follows: 

• The statistical functions n^iq) and 1 + n^iq) are replaced by their classical limits T/q; 

• Dispersion relations and propagators use pi, qi in place of pj, qi. 

Furthermore, at the order considered, the diagram is dominated by small exchange momen- 
tum and we can approximate the argument of the 5- function as {p — q)"^ —if = 2p-q+0(p?), 
where <?« = ^ sin(agj). Thereby the lattice version of Eq. (|3.6p becomes 

ATrg^T^CpNc d^p f^'" d^q 5(p • q) 

'^latt — 



3 (27r)3 (p2 + ^2^^^^^J2 (2^ 



^3 



3vr J_^/a (2vr)3 ^^2^2 \^ amo.iatt 

Here we made use of the fact that for p ~ Tn^ iatt ^ l/oi P can be viewed as a continuum 
variable, so one can carry out an angular integral to remove the (^-function. We calculate the 
constant accompanying the logarithm, denoted by ... in Eq. (j3.7p . in Appendix A, finding it 
to be 1.8313(2). 
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In Eqs. (|2.ip . (|2.2p . we saw the correspondence 



where the defining expression for T, reads 

S _ r/"" d^q 1 



ATra - J_^/, (2vr)3 g2 ' (3-9) 

Comparing the coefficients of the logarithms in Eqs. ()3.6p . (j3.7p . on the other hand, suggests 
the correspondence 

, . (3.10) 

The difference between Eqs. ()3.8p . ()3.10p is a manifestation of the ambiguity in the matching 
of the continuum and lattice theories that was mentioned in Sec. El More generally, noting 
that the dispersion relation in Eq. (13. 5p gives a group velocity 




and defining 



we can write 



(3.11) 

_ d3q \v^ronMT, (3 12) 



4™ (2^)3 g2 



Kiatt ~ ^ ^ X In h 1.831 . (3.13) 

A similar correspondence was found in ref. [28j for a number of other quantities: scaling away 
Yatf Debye screening length involves Yi^o = S: infrared magnetic damping involves S„i ; 
and the plasmon oscillation frequency involves o What we have shown is that Coulombic 
scattering (electric damping) involves S^-i. The numerical value of each S^n is given in 
Table [TJ The fact that they do not coincide means that there is some ambiguity in how to 
relate the lattice and continuum theories. Nevertheless, the general structures of the answers 
in Eqs. (13. 6p . ()3.7p . including the existence of logarithms, are the same, whereby we can 
conclude that the dynamics of the two theories indeed bear a strong qualitative resemblance 
to each other. 



^Each V dependence arises from simple physics. Kjatt mean squared momentum a charged particle 

absorbs due to Coulomb interactions with passing excitations. An excitation with velocity v has a flux factor 
suppressed by v but it interacts for 1/v times as long, giving a force-squared enhanced by v~^; hence Kiatt 
involves . Magnetic damping is similar but it involves magnetic forces, which are suppressed by v relative 
to Coulombic forces. Therefore the mean squared momentum exchange scales as v''\ leaving only the flux 
factor . Plasma oscillations involve the mean squared current generated by an oscillating electric fleld; the 
current is proportional to v of the charges, so ojpj ~ . Debye screening is a thermodynamic property so it 
shows no scaling with group velocity. 
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"standard" [14] 


"improved" [33] 




1.6222746498 


1.78576519 




2.1498783949 


2.13792379 




3.1759115356 


2.783189232 




5.5079614967 


4.1679252 



Table 1: Values of for the standard and "improved" lattice actions. The spread of values is a 
measure of how different the structure of the lattice Hard Thermal Loops is from the continuum ones. 



The discussion so far has assumed that we are in the weak-coupling regime, i.e. that /3l 
is large. If /?l decreases so much that all three momentum scales are of the same order of 
magnitude, then the behavior of Kiatt changes significantly. Scaling Kiatt dimensionless by 
multiplying with a^, and making use of the definition of the weak-coupling behavior in 
Eq. p.l3p corresponds to 

a^^/^iatt 4ln/3L , (3.14) 



while at small /3l we find, through the arguments in Appendix B, the behavior 

3 '3l<1 1 , . 

a Ki^tt ~ ■ (3.15) 

The physics behind this functional form is very specific to the nature of the lattice variables, 
however, so we do not expect any analogy with the continuum quantum theory in the latter 
regime, and refrain from a further discussion here. 



4. Numerical Results 

In the previous section we have verified that, at weak coupling, the electric field correlator 
in classical lattice gauge theory behaves quite similarly to physical QCD: if we fix the lattice 
spacing a by equating the Debye screening lengths, the leading-logarithmic k of the classical 
lattice theory is larger than that in the quantum continuum theory by a factor S^-i /S^o ~ 
5/3. We now proceed to larger values of the coupling (smaller /3l) with the help of numerical 
simulations. At relatively small coupling (large /3l), we can check how fast the weak-coupling 
regime is approached. At stronger coupling (intermediate /3l ~ 1), we can find out whether 
the leading-order weak-coupling result is reasonable even in order of magnitude, and whether 
it underestimates or overestimates the actual behavior. This then gives us some guidance for 
what to expect in QCD. 
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SU(2) 



SU(3) 




Figure 2: Numerical results for the intercept Ki^tt (open symbols), compared with the weak-coupling 
prediction from Eq. ([XT^ (line). The left plot is for SU(2), the right one for SU(3). Note that 
1//?L = g^Ta/2Nc scales like as, assuming the matching a ^ 1/T (cf. discussion after Eq. (|2.2p ). 



4.1. Intercept at a; ^ 

Our numerical results for Kiatt) compared with the leading-order weak-coupling result, are 
shown in Fig. [21 both for SU(2) [included because a large /3L-range could be scanned with a 
modest numerical effort] and for SU(3)S We note, first of all, that at large the results 
approach the analytic ones of Eq. (I3.13p . However, as soon as /3l ^ 100, the non-perturbative 
results deviate from the leading-order ones. The non-perturbative results are always larger 
than the perturbative estimate. For /3l = 1...10, a crossover takes placed from one type 
of behavior to another. At /3l ^ 1) the results approach the behavior of Eq. (|3.15p . (We 
have not worked out the numerical prefactor for Eq. ()3.15p . and hence do not show the 
corresponding curves in Fig. O) 

In order to make quantitative use of the numerical results, it is convenient to change the 
units of both axes. Recalling the definition of j^^^.^. from Eq. (j2.2p . we choose the variable 
g^NcT/rriD ig^tt = 2Nc{ii /T.Pi^Y/'^ as the x-coordinate; this quantity is the ratio of the g^T 
to gT scales and is therefore the expansion parameter for perturbation theory at the scale 
gT. We also divide a^^iatt by the coefficient of the leading logarithm, a^g'^C-pTrri^Y^^^/Q-K = 
CfA^^E/Stt^/Sl- In these units, the lattice results have a direct counterpart in the continuum 

''in the numerical implementation the theory is discretized in time as well as space, but with a much finer 
spacing, and our numerical results for Kj^tt represent the limit of zero temporal spacing. We have also checked 
that our results contain no significant finite volume or non-zero ua artifacts. 

^We have checked that there is no actual phase transition in the thermodynamics of the system. 
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theory. The weak-couphng regime is plotted in the new units in Fig. [3l 

In the continuum theory, corrections of 0{g^) to Kcont have recently been determined [8], 
and it is now interesting to compare the results. According to ref. [8], 

In h Ccont + -Dcont V ■■■) , (4.1) 



^cont 



67r V mD,cont "lD,cont 



with -Dcont = 0.7767. The physics giving rise to -Dcont involves only the length scale gT 
and should be reproduced on the lattice; however, it depends on the structure of the Hard 
Thermal Loops in an essential way, so the lattice value could differ by up to 0(50%) as 
discussed in the previous section. Still, this motivates a fit of the lattice data to the form 

Klatt = ^ ■ — ^ In h Clatt + ^latt h • • • , (4.2) 

OVr \ h amD,latt "iD.latt / 

where Ciatt = 1-831 x S^-i/S = 3.176. The coefficient Aatt can confirm the sign and 
approximate magnitude of -Dcont- It can also tell us about the next terms in the expansion, in 
particular whether the 0{g^) calculation is an underestimate or an overestimate of the real 
K. Note, however, that the next term, of 0{g^), would receive contributions not only from 
the scale gT but also from the scale vrT, so the relation between the lattice and continuum 
theories becomes less precise at this order. 

Fig. [3] shows the (1-parameter) fit to the lattice data according to Eq. (j4.2p . The fit is very 
good out to g'^NcT/mD^ig.tt ~ 1-5. We extract the value -Diatt = 0.87(4), in surprisingly good 
agreement with Dcont = 0.7767; and we see that the same coefficient -Diatt fits the SU(2) and 
SU(3) data, just as the continuum computation predicts. However, at larger couplings Kiatt 
rises above the fitted behavior, particularly for the group SU(3). 

The above results no doubt depend on the details of the numerical implementation of the 
lattice theory and of the electric field operator. As a check on the robustness of our results, 
we re-compute them using the "improved" lattice action of ref. [33]. This action is tree-level 
improved so that the IR behavior naively coincides more tightly with the continuum, as shown 
for instance by the better IR behavior of the free-theory correlator k{uj) in Fig. [TJ However 
the UV behavior still has (different) anisotropic non-ultrarelativistic dispersion, so the Hard 
Thermal Loop effects are not those of the continuum (though they are somewhat closer, as 
reflected by the slightly narrower spread of the E^n values in Tabled]). Therefore it is better 
to think of this implementation as "different" rather than truly "improved." 

Fig. Ulshows k{0) as a function of the lattice coupling for the "different" / "improved" lattice 
action for the group SU(2). While the lattice constants Ciatt = 2.5 and -Diatt = 0.64 (this 
time both are fitted) differ from the "standard" action values, the qualitative message is the 
same; at weak coupling the behavior appears to be well described by next-to-leading order 
perturbation theory, but at stronger coupling perturbation theory is an underestimate. But 
while for /3l ^ 1 the two implementations give similar qualitative results, at extremely strong 
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g^NT/m^, „ 

° c D.latt 

Figure 3: Ki^tt, normalized to the leading-order perturbative behavior, expressed as a function of 
the expansion parameter 5^-/VcT/TOD,iatt related to corrections from the scale gT. At weak coupling, 
the SU(2) and SU(3) results agree and are well fit by the 0{g^) perturbative behavior. At stronger 
coupling, Kiatt rises above the perturbative fit, by a group dependent amount. Compared with Fig. [2l 
the horizontal axis is restricted to < 0.77 for SU(2), < 0.34 for SU(3). 



coupling the behaviors are not even qualitatively the same. This reinforces our belief that 
the /?L ^ 1 behavior is a lattice artifact with no bearing on QCD. 

We finally attempt a rough order-of-magnitude estimate for which /^L-range corresponds 
to the situation met in heavy ion collision experiments. Combining the matching from below 
Eq. (12. 2p with the definition of we get 

2N, + Nf 

77^ • 4.3) 

The relevant value of Os can in the present context probably best be approximated by taking it 
from the dimensionally reduced effective theory [M] , to which the classical lattice gauge theory 
reduces in the case of equal-time observables. In this limit the coupling has been computed 
up to 2-loop level [30]; for A^f = 3 the values are ~ 3. ..2 for T/A^jg 1...4, corresponding 
to tts = 0.24. ..0.16, and subsequently /3l = 2. ..3. This corresponds to g^A^cT'/w-D.iatt = 4.. .3. 

Remarkably, in the range /3l = 2... 3, the numerical SU(3) values in Fig. [2] exceed the weak- 
coupling result by as much as an order of magnitude! Though these values of /3l are so small 
that the matching cannot be trusted on any kind of quantitative level, such a huge effect is 
still encouraging both from the experimental point of view [35], where the apparently very 
rapid thermalization of heavy quarks remains a mystery, as well as from the point of view of 
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SU(2) SU(2) 




Figure 4: ^i^tt using the "improved" lattice action. Left: the overaU behavior in lattice units. Right: 
a magnification of the weak-coupling regime, normalized to the leading-order perturbative behavior. 
The weak-coupling behavior is in good qualitative accord with the standard action, but the strong 
couphng behavior is qualitatively different (cf. Figs. [H [3]). 



following the suggestion of ref. [22] in order to measure Atcont with Euclidean lattice Monte 
Carlo methods. Indeed, there may well be an exciting qualitative discovery to be made on 
the lattice. 



4.2. General shape of the spectral function 

On the point of lattice Monte Carlo simulations, ref. [22] argued that the Euclidean analogue 
of Eq. (jl.ip leads to a correlator, denoted by Ge{t), which has a non-trivial continuum limit 
and can be related to the intercept Kcont through standard relations. Specifically, the task 
would be to invert the relation 

roo J cosh — t) LO 

GE{r) = / — K,ont(o.) ^ , ' . (4.4) 



- cosh^ 

It is a problem, though, that strictly speaking the relation in Eq. (j4.4p is not invertible 
without further input. In practice, this means that a certain Ansatz (sometimes called a 
prior) is needed, which is then refined through the numerical data. For this reason, significant 
efforts have been devoted to analytic computations of spectral functions in the presence of a 
spatial lattice, in the limit of a high temperature, for cases such as the 2-point correlator of 
the vector current of heavy quarks [36] . 

We can now use our data, both perturbative as well as non-perturbative, to obtain an 
Ansatz for the spectral function Kcont('^)- In Fig. jSj results are shown for the function 
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Figure 5: Numerical results for the function a^Ki^tt (open symbols), compared with the weak-coupling 
prediction from Eq. for /3l = 24 (line). The left plot is for SU(2), the right one for SU(3). For the 
(unphysical) case /3l = 1, included as a reference for the discussion in Appendix B, we have divided 
the central values (but not the error bars) by a factor 10. 



a^Kiattlf-^) at various together with a comparison with the free theory result. Noting that 
on the 4-dimensional lattice, (3 = Nrar, where A'^,-, ar are the number of lattice points and 
the lattice spacing in the time direction, respectively, and naively enforcing the replacement 
of the classical limit of the Bose-Einstein distribution function, T/u, by the corresponding 
quantum mechanical expression, 1/2 + 713(0;), we can expect Kcont(w) to behave as 

LOarNr ^. ( LOarNr\ , . , . 

'*cont(w) ^ 7: COth Ulatt('^) ■ (4.5) 



In particular, for uj <^ T, Kcont('^) should be completely flat just like Kiatt('^); moreover, in 
general, Kcont('^) should show no peaks other than at a; ~ (1.5 — 3.0)/a, where a is the spatial 
lattice spacing. We consider these qualitative features to be relatively robust, and they can 
in any case serve as crosschecks on particular practical inversions of Eq. ()4.4p . 

Finally, we remark that the corresponding spectral functions computed for J\f = 4 Super- 
Yang-Mills theory at infinite 't Hooft coupling in continuum show an analogous smooth 
behavior at small frequencies, taken over by ultraviolet physics at u; ~ T [31 137j . 



5. Summary and Outlook 

The purpose of this paper has been to make use of classical lattice gauge theory, in order to 
gain insights on the dynamics of QCD in the temperature range accessible to current and near- 
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future heavy ion collision experiments. We have stressed, in particular, that classical lattice 
gauge theory is a multiscale system just like QCD; unlike QCD, however, it easily lends 
itself to non-perturbative simulations of real-time observables, in both the weak-coupling 
and strong-coupling regimes. Thereby a semi-analytic understanding can be obtained of 
many interesting observables, without changing the number of color degrees of freedom or 
introducing unphysical infrared fields. 

More specifically, we have elaborated on the heavy quark momentum diffusion coefficient, 
denoted by k, which determines the heavy quark thermalization rate through linear response 
relations. This quantity belongs to the general class of observables which are "dominantly" 
influenced by momenta around the Debye scale, p ~ gT. We have shown explicitly through 
a weak-coupling analysis that while the physics of the classical lattice gauge theory differs 
from that in QCD on the quantitative level, by effects of up to 50%, the qualitative features 
of the dynamics do remain intact. 

Proceeding from weak coupling towards intermediate coupling, we have furthermore shown 
that the leading-order weak-coupling expression, and even the larger next-to-leading order 
expression, underestimate the non-perturbative result. Given the close analogy with QCD, 
the same statement should be true on that side. This seems to give realistic hopes that a 
future quantitative determination of k through 4-dimensional lattice Monte Carlo simulations 
will reveal a large thermalization rate, which might help to explain the surprisingly rapid 
thermalization that has been observed at RHIC experiments [35j . 

Finally, with regard to Monte Carlo simulations, we have explored the general structure 
of the spectral function corresponding to the Euclidean electric field correlator that can be 
used for determining the thermalization rate [22]. We find that apart from a single peak at 
the scale of the spatial lattice spacing, the spectral function has little structure, both at weak 
and at intermediate coupling. This should be an encouraging message with respect to the 
analytic continuation needed in the analysis of the Monte Carlo simulations. 
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Figure 6: The diagrams contributing to /ti^tt in Coulomb gauge. The straight hnes represent temporal 
Wilson lines, the closed squares electric fields, and the grey bubble the gauge field self-energy. 

Appendix A. Weak-coupling regime in classical lattice gauge theory 

We give in this appendix some details concerning the computation of the constant accompa- 
nying the logarithm in the leading-order weak-coupling result, Eq. (j3.7|) (or Eq. (j3.13p ). 

In the "standard" implementation [T3] of classical lattice gauge theory, where time is con- 
tinuous and Minkowskian, the electric field strength has the form 

aE.i.) ^ -l[W.(.)lCt(.) + M.) - U.i.)M. + a}>U} (.) . (A.1) 

Here Ui are the spatial link matrices. For a perturbative computation we write Ui = 
exp{iagA^T''), where are Hermitean and assumed normalized as Tr [T^T'^] = As 
usual [38], Fourier representations of the spatial variables are most conveniently chosen as 
A^{x) = Jj^ A^{K)e^^'^^~^"'^^'^h Furthermore a gauge needs to be fixed; like in the continuum 
computation [2], it is convenient to choose a Coulomb gauge so the propagator splits into a 
transverse spatial part and an Aq field propagator which has no on-shell spectral weight. In 
this gauge the result emerges from two graphs, depicted in Fig. E) the self-energy diagram, 
mediated by the vertex 

A^oiP) A'^iQ) A'jiR) 5{P + Q + i?)|r'"^5i,- cos(aK/2)(i?o - Qo) , (A.2) 

plus from one additional graph, namely the bubble diagram sourced by the second term in 
the expansion of the first term of Eq. ()A.ip . 

aE,{x) = - + \ a'gT'f''''[dtA1{x)]Af{x) + ... . (A.3) 
The final result can be written as 



'^latt 



3 (2^)3 (p2+ ^2^)2 y_/j2vr) 



\3 



where the limit mo *C 1/a is assumed. 



„2 JL ^pjqf +pj{pi - qif p^qfipi - qif^ 



cp + q^ 



(A.4) 
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For the numerical evaluation of Eq. ()A.4p , one can for instance integrate explicitly over one 
of the momentum components, to remove the (^-function, and carry out the remaining five- 
dimensional integral numerically, simplifying the range by making use of various symmetries. 
More refined strategies are certainly possible but not necessary if only a few digits are needed. 



Appendix B. Strong-coupling regime in classical lattice gauge theory 

Figure [2] shows that, in the limit /3l — > 0, the electric field autocorrelator diverges as fi-J'^'^-, 
while Fig. [5] shows that the frequency spectrum for Kiatt(i^) becomes tightly peaked at small 
frequencies. What is going on in this regime, and could it have anything to do with QCD? 
Here we show that the answer to the latter question is almost certainly negative. 

To do so we need to discuss a few features of the numerical simulation, which we have 
otherwise left to the references. Fixing to the temporal gauge, which is convenient because 
the temporal Wilson lines in the definition, Eq. (jl.ip . are identity operators, the continuum 
Yang- Mills theory is described by gauge fields j4j(x, t) and their canonical momenta, the 
electric fields Ei{x,t). On the lattice, the degrees of freedom are the dimensionless electric 
fields £i = a'^gEi+0{o?) and the gauge links, Ui{x) = eyiY){iagAi{x)). We can write £i = S^T'^, 
where are Hermitean generators of the SU(A^c) algebra, normalized as Tr [T^'T'^] = b^'^ jl. 

The lattice simulation proceeds by sampling initial configurations with a classical Hamilto- 
nian, and then evolving the configurations in real time through classical equations of motion. 
The time evolution of the link matrices is 

adtV,{x)=iE,{x)Vi{x) , (B.l) 

while for the electric fields it is 

adtS-ix) = 2^ImTr|r^ Uj{x)Ui{x + aj)U^j {x + ai)U} {x) 



+ U]{x - aj)U,{x - aj)Uj{x + ai - aj)ul{x)\^ } . (B.2) 



What is important here is that the spatial link variables are compact. Therefore the size 
of the time derivative of is bounded. At sufficiently small /3l (strong coupling) the link 
matrices become essentially random elements of the group, and the typical size of \adt£i\ 
saturates. On the other hand, the mean-squared value, P) ~ 1//?L) does not saturate but 
increases linearly as is made large. 

In this regime, the electric fields feel an essentially random force of fixed mean-squared 
value, and evolve much like heavy particle velocities in classical Langevin dynamics. The time 

, 1 /2 

scale for the link matrices Ui to rotate by an 0(1) angle is t/a ~ l/|<?f | ~ (3-^ ■ The force on 
the electric field, Eq. (|B.2p . involves a product of four links which each rotate independently 
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at (generically) irrationally related frequencies; the product of four such randomly rotating 
group elements should show no periodicity or quasi-periodicity. Therefore the coherence time 
of the random force on £■' is set by the time for a link matrix to rotate by an 0{1) angle, i.e. 
just t/a Ph"^ ■ A random variable £^ with mean squared value \£i\'^ ~ experiencing 
a random force of magnitude ~ 1 with a coherence time ~ /J^^^ behaves as {£^{t)£j'[0)) 

exp(— |t|/T) where r ~ Integrating over t and inserting (|<?fp) ~ P^^, we 

—5/2 

conclude that k ~ /J^ for /3l <C 1. This description also predicts that the support of 
^iatt(i^) should become narrow with width aw ~ a/r ~ /3-^ . 

However we emphasize that this behavior is an artifact of the electric fields being non- 
compact while the gauge links Ui are compact. Such a disparity is absent in the quantum 
theory so the effect is an artifact of the classical lattice discretization. Concretely, we find 
a different qualitative behavior in the small-/3L regime of the "improved" description (cf. 
Fig- SI)- Therefore we believe that the behavior of Ki^tt in tbe small- /?l regime has nothing to 
do with real QCD. 
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